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Many-body quantum coherence and interaction blockade in 
Josephson-linked Bose-Einstein condensates 
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Abstract. - We study many-body quantum coherence and interaction blockade in two Josephson- 
linked Bose-Einstein condensates. We introduce universal operators for characterizing many-body 
coherence without limitations on the system symmetry and total particle number N. We reproduce 
the results for both coherence fluctuations and number squeezing in symmetric systems of large N, 
and reveal several peculiar phenomena that may occur in asymmetric systems and systems of small 
TV. For asymmetric systems, we show that, due to an interplay between asymmetry and inter- 
particle interaction, the coherence fluctuations are suppressed dramatically when \Ec/Ej\ <C 1, 
and both resonant tunneling and interaction blockade take place for large values of \Ec/Ej\, 
where Ec and E.j are the interaction and tunneling energies, respectively. We emphasize that the 
resonant tunneling and interaction blockade may allow creating single-atom devices with promising 
technology applications. We demonstrate that for the systems at finite temperatures the formation 
of self-trapped states causes an anomalous behavior. 



Utilizing the current experimental techniques, it is 
possible to prepare atomic Bose-Einstein condensates in 
which all atoms occupy the same single-particle state. To 
explore the quantum coherence between two such conden- 
sates, it is natural to link them through the Josephson 
tunneling [1], and then release for producing interference. 
Atomic interferometry [2] provides an important experi- 
mental tool for the fundamental studies in atomic, optical, 
and quantum physics, as well as promising applications 
in high-precision measurements and quantum information 
processing [3]. By loading condensates into double- well 
potentials, atomic interferometers have been realized ex- 
perimentally [4-7]. Furthermore, some many-body quan- 
tum effects such as coherence fluctuations [6,8], condi- 
tional tunneling [9], number squeezing [10], and long-time 
coherence assisted by number squeezing [11] as well as 
finite-temperature effects [8] have been observed in exper- 
iment. However, clear theoretical explanation for many of 
those effects, such as the full physical picture of the coher- 
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ence fluctuations [6,8], the resonant tunneling and the in- 
teraction blockade [12-15] are still largely missing. Below, 
we present a general approach for analyzing many-body 
quantum coherence and interaction blockade and give a 
foil picture of the coherence fluctuations. Particularly, our 
results of resonant tunnelling can qualitatively explain the 
conditional tunneling which has been observed recently 
[9]. 

To study quantum fluctuations in the Josephson-linked 
coherent systems such as Bose-Einstein condensates in 
double- well potentials, one has to analyze the system 
within the full quantum theory. Under the conditions of 
tight binding, the quantum system can be mapped onto 
a two-site Bose- Hubbard model [10]. Then, the num- 
ber fluctuations can be easily analyzed by applying the 
techniques well- developed in quantum optics [16]. Where, 
however, the coherence fluctuations could not be clearly 
depicted. Previously, the phase fluctuations have been 
used as a measure characterizing the coherence fluctua- 
tions. However, due to the absence of the phase operators 
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in the Bose-Hubbard Hamiltonians, one has to introduce 
the quantum-phase models [17,18]. The quantum-phase 
models describe qualitatively the systems with small num- 
ber differences, but they can not describe the systems of 
large number differences and fail to maintain the hermitic- 
ity for asymmetric systems. Additionally, because the 
phase is measured modulo 2ir, the fluctuations of phase 
do not provide the best quantitative characterization for 
the coherence fluctuations [11]. 

In this Letter, we introduce Hermitian operators for the 
phase coherence rather than a relative phase, and then 
analyze the coherence fluctuations and interaction block- 
ade. The coherence operators are defined with the atomic 
creation, annihilation, and number operators, so that they 
can be calculated directly with the general Bose-Hubbard 
models. The operational approach has two advantages: 
first, its validity is independent of the total particle num- 
bers and, second, it works for both symmetric and asym- 
metric systems. We examine how the coherence fluctua- 
tions depend on the total number of particles, the inter- 
particle interaction energy, the junction tunneling energy, 
asymmetry, and temperature. Without loss of general- 
ity, for all calculations about the fluctuations in zero- 
temperature systems, we only consider the ground states 
of the full quantum model: quantized Bose-Josephson 
junction. For the calculations about the fluctuations in 
finite-temperature systems, we use both the full quantum 
and the mean-field models. 

Below, we not only reproduce the pictures of coher- 
ence fluctuations and number squeezing in large-size sym- 
metric systems which can be described by the quantum 
phase models, but also reveal several peculiar phenomena 
in small-size systems and asymmetric systems which can 
not be described by the standard quantum phase mod- 
els. Such as, (i) due to an interplay between the system 
asymmetry and inter-particle interaction, single-atom res- 
onance effects suppress the quantum fluctuations of co- 
herence and the interaction blockade occurs between the 
neighboring resonant peaks; and (ii) the occupation of 
highly-excited states of degeneracy (macroscopic quantum 
self-trapping) causes an anomalous behavior of the coher- 
ence fluctuations at finite temperatures. These results are 
qualitatively consistent with the recent experiment data 
on many-body quantum phenomena [6,8,10,11] and finite- 
temperature dynamics [8]. In particular, the resonant ef- 
fects are confirmed by the most recent experimental ob- 
servations of the conditional tunneling [9]; the physics of 
the interaction blockade and resonant tunneling open a 
route for designing single-atom devices for applications in 
sensitive metrology and information technology. 

We consider the quantum model of two Joscphson- 
linkcd atomic condensates with the conserved total num- 
ber of particles. From the physical point of view, this 
situation can be realized by confining atomic condensates 
in a double- well potential [4-7]. Eliminating the constant 
terms including the total number of atoms, the system 



Hamiltonian 

5 E 
H = -J(a^a 2 + a 1 aj) + -(n 2 -n 1 ) + —(n 2 -n 1 ) 2 (1) 

describes a quantized Bose-Josephson junction. Here, 
a~j(aj) are the atomic creation (annihilation) operators 
for the j-th (j = 1,2) well, and rij are the corresponding 
number operators. Thus, the system behavior is deter- 
mined by the total number N, the charging (inter-particle 
interaction) energy Eq, the inter- well tunneling strength 
J, and the asymmetry parameter S. It has been suggested 
that the symmetric systems with negative charging ener- 
gies can be used to realize the Heisenberg-limited inter- 
ferometry [19]. Here, we consider both symmetric and 
asymmetric systems with the positive charging energies. 

Under the mean-field approximation, aj — ► y/nj e l< ^ j 
and — > ^/nj e~' l ^ j , the macroscopic quantum behav- 
ior (such as Josephson oscillations, self-trapping [20], and 
dynamical bifurcation [21]) obeys a classical non-rigid pen- 
dulum Hamiltonian, 

H cl = E c n 2 /2 + 5n- Ej\Jl - (2n/N) 2 cos <j), (2) 

with the quasi-angular momentum n = (n 2 —n\)/2 and the 
relative phase <f> = (fi 2 — 4>i- Here, Ej — NJ denotes the 
junction energy. For the symmetric systems (S = 0) with 
small number differences (\2n/N\ <~ 0), one can introduce 
the quantum phase model via n — ► i-j^ and (2n/N) — ► 
(see Ref. [17,18] and references therein). However, it can 
not depict the behavior of systems with large number dif- 
ferences or self-trapping, and it fails to maintain hcrmitic- 
ity when 5^0. 

To explore the coherence fluctuations, we introduce op- 
erators for the phase coherence. Drawing upon the quan- 
tum phase concept [22] for single-mode coherence, we de- 
fine the operators for two-mode first-order coherence as 



(a 2 ai + a 2 af) 



COS0 = 

sin0 = 



y/2 (2 ni n 2 - 


-ni + n 2 ) 


{(a^ai — 


a 2 af) 



(3) 



^2 (2nin 2 + m + n 2 ) 



These two operators satisfy the conditions 
(sin 2 4> + cos 2 4>) = 1- Without considerations in to- 
tal particle number conservation and inter-particle 
interaction, similar operators have been used to study the 
phase coherence of photons from different lasers [23] . So 
that the previous results of photons can not be applied 
to systems of conserved total particle numbers and 
inter-particle interaction, such as the Josephson-linkcd 
condensates considered in this paper. Similarly, one 
can define the operators for two-mode second-order 
coherence, cos(2(/>) = [(a^ai) 2 + (a 2 af ) 2 ] /Kg and 
sin(2(/>) = i [(a^ai) 2 — (a 2 af) 2 ] /Ks, where the constant 
Ks = (2nin 2 + n\ +n 2 ). Fluctuations of the two-mode 
first-order coherence arc quantitatively characterized by 
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the variance 

A(cos0) = (cos 2 <f)) — (coscj) 
Here, the expectation for cos <fi is 



(cos <f>) 



a/2 (2nin 2 + n\ + n 2 ) 
and the expectation for cos 2 </> is 



(cos 2 (f>) 



1 | ((g 2 + fll ) 2 ) + ((a 2fl +) 2 ) 

2 2 (2nin 2 + rii + n 2 ) 



(4) 



(5) 



(G) 



The variance A (cos <fi) vanishes when two condensates be- 
comes phase-locked. Otherwise, the variance approaches 
the value 1 /2 when two condensates have no phase correla- 
tions. The above variance has been successfully employed 
for explaining the q- vortex melting [24]. 
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Fig. 1: Coherence fluctuations and number squeezing in sym- 
metric Bose-Josephson junctions, (a) Coherence fluctuations 
A(cos<j!>) vs. Ec/Ej; and (b) inverse square of number squeez- 
ing 1/S 2 vs. N 2 Ec/Ej. Dashed lines, dashed-dot lines, cir- 
cles, and dots correspond to the cases N = 2,4, 50 and 100, 
respectively. For sufficiently large N, the curves approaches 
the result obtained from the Gaussian ansatz applied to the 
probability amplitudes (solid line). 

For symmetric systems at zero temperature, the co- 
herence fluctuations A(cos</>) and the number squeezing 
S = y/N/AN r (where AN r = (N?) ~ (N r ) 2 is the vari- 
ance for relative numbers N r — n 2 — ni) show differ- 
ent parametric dependencies. The coherence fluctuations 
grow monotonously with the ratio Ec/Ej = Ec/(NJ), 
however, the number squeezing grows monotonously with 
the ratio N 2 E C /E,j = NE C /J. The behavior of the co- 
herence fluctuations is qualitatively consistent with the 
experimental observations made via adjusting the sepa- 
ration distance [6] or the barrier height [8], which corre- 
sponds to varying Ej. Because the coherence time grows 



with the number squeezing [25], the different paramet- 
ric dependence of the coherence fluctuations and num- 
ber squeezing indicates that one can increase the coher- 
ence time without loss of coherence, by increasing the 
total number N or with some loss of coherence, by in- 
creasing the ratio Ec/Ej. This result justifies qualita- 
tively the experimental approach for adjusting the coher- 
ence time via varying the barrier height (i.e., controlling 
Ej) [11], in which the coherence time grows with the bar- 
rier height. For sufficiently large N , the parametric depen- 
dence of coherence fluctuations and number squeezing are 
consistent with the previous results obtained from other 
approaches [f7,25,26]. Our operational approach also ex- 
plores the finite-size effects in systems of small N. For 
commensurate systems (whose N are even integers) of very 
weak or strong interactions, the coherence fluctuations are 
almost independent of N . However, in the intermediate 
region of Ec/Ej ~ 1, the coherence fluctuations decrease 
with N. In Fig. I, we show the parametric dependence 
for the ground-state coherence fluctuations and number 
squeezing for symmetric systems with different total num- 
ber of particles. 

For asymmetric systems at zero temperature and with- 
out inter-particle interaction [Eq — 0), both coherence 
fluctuations and number squeezing increase monotonically 
with the asymmetry parameter S. For small values of 
Ec/Ej, before the asymmetry S starts dominating the 
dynamics, the coherence fluctuations can be suppressed 
by an interplay between asymmetry and repulsive inter- 
particle interaction; the suppression becomes more sig- 
nificant for larger values of Eq- This suppression of the 
coherence fluctuations originates from the different local- 
ization mechanisms dominated by asymmetry and repul- 
sive inter-particle interaction. The asymmetry tends to 
localize all atoms in the lower well, however, the repulsive 
inter-particle interaction (positive charging energy) tends 
to localize atoms in two wells without number difference. 
In Fig. 2, we show the coherence fluctuations and number 
squeezing in asymmetric systems with N = 100, Ej — 100 
and Ec/Ej < 1. 

For asymmetric systems with large values of Ec/Ej 
at zero temperature, the coherence fluctuations can be 
suppressed resonantly, due to degeneracy and quasi- 
degeneracy of the energy levels occurring at some partic- 
ular values of the ratio S/Ec- For systems without tun- 
neling, the eigenstate of quasi-angular momentum n has 
the energy E*(n) = 5n + Ecn 2 /2, and so that two neigh- 
boring levels become degenerate when E*(n) = E*(n + 1). 
This means that, near the values of 8/Ec = — (n + 1/2) 
for n = (-N/2,-N/2 + l,---,N/2 - 1), at least one 
atom can tunnel freely between two wells, even in the 
limit of very weak tunneling (i.e., for very small values 
of Ej). Therefore, the coherence fluctuations from inter- 
particle interaction will be suppressed resonantly, with 
the resonance width depending on Ej. This is a single- 
atom resonance effect caused by the degeneracy and quasi- 
degeneracy between the states of the neighboring angular 
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Fig. 2: Coherence fluctuations and number squeezing in asym- 
metric systems with TV = 100, Ej = 100, and E c /Ej < 1. 
(a) Coherence fluctuations vs. S/Ej; and (b) inverse square of 
number squeezing 1/S 2 vs. S/Ej. 

momenta n. For very small values of Ej/Ec, we perform 
a perturbation analysis by assuming that the perturbed 
ground state is a superposition of the two lowest states of 
the system without tunneling. Namely, the ground state 
\"4>)q = ot\ \L\) + a 2 \L2)i where \L\) and \L 2 ) are the two 
lowest states. The two coefficients are determined by the 
lowest energy state of the following eigenvalue equation, 

( H u flia \ ( ax \ „( «i \ 
^ H21 H 22 )\a 2 ) \ a 2 J ' 

in which Hij = (Li\H\Lj) = 1,2). The perturba- 
tion results confirm the resonant suppression of coherence 
fluctuations, and they show that the resonance width in- 
creases with Ej/Ec- Similarly, if we prepare a state of 
a quasi-angular momentum n, the correlated tunneling of 
an atomic pair (two-atom resonance) should correspond to 
the degeneracy condition E*(n) = E*(n + 2). Our results 
for single- and two-atom resonances explain consistently 
the recent observations of conditional tunneling [9]. In 
Fig. 3, we show the coherence fluctuations in the ground 
states of asymmetric systems with N = 100, Ej = 100 and 
large values of Ec- For a very large value of the charging 
energy, Ec — 800 (green lines), our numerical results differ 
slightly from the the results calculated by the perturba- 
tion analysis. Even for a moderate value of the charging 
energy, Ec — 100, our numerical results clearly show the 
resonant suppression of the coherence fluctuations. Ap- 
parently, the resonance for the case Ec = 100 is wider 
then that for the case Ec = 800. 

The single-atom resonant tunneling is a signature of in- 
teraction blockade induced by the s-wave scattering be- 
tween atoms. It is well known that the Coulomb blockade 
of single-electron tunneling occurs because of the charg- 
ing energy [12]. In the Bose-Josephson junction (1), the 
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Fig. 3: Coherence fluctuations in the ground states of asym- 
metric systems with TV" = 100, Ej — 100 and large values 
of Ec- (a) Coherence fluctuations A(cos</)) vs. 5/Ec; and 
(b) the enlarged region between two dot lines shown in (a). 
Blue and green lines correspond to the values Ec = 100 and 
Ec = 800, respectively. Resonant suppression of the co- 
herence fluctuations occurs at S = —Ec(k + 1/2) for k = 
(-N/2,-N/2 + 1, ...,N/2 - 1), marked by dashed red lines 
in (b). 



asymmetry and the inter-particle interaction act the roles 
of the applied voltage and the charging interaction, re- 
spectively. As a results of the interaction blockade, adi- 
abatically sweeping the asymmetry, atoms in the ground 
state tunnel through the potential barrier one-by-one and 
the coherence vs. asymmetry dependence is a set of sharp 
peaks between plateaus of interaction block. The interplay 
of the single-atom resonant tunneling and the interaction 
blockade allows one to create so called single-atom de- 
vices, which would provide promising application in sensi- 
tive metrology and information technology like their elec- 
tronic counterparts [14, 15]. For a strongly interacting 
system in optical lattices with asymmetry imposed by 
a harmonic potential, the appearance of Mott-insulator 
shells [27, 28] is also the results of interaction blockade. In 
Fig. 4, we show resonant tunneling peaks and interaction 
blockade plateaus for a quantized Bose-Josephson junction 
of N = 10, Ej = 10 and E c = (1000, 100, 10). In which, 
A(7V r ) and < N r > are the quantum fluctuation and ex- 
pectation of the relative particle number N r , respectively. 
The regular steps in < N r > and sharp peaks in A(N r ) 
clearly show the single-atom resonant tunneling and the 
interaction blockade. It also show that the interaction 
blockade becomes more and more significant when the in- 
teraction energy Ec increases. In the recent experiment 
of atoms confined in double- well lattices [9], N = 2 and 
Ej / Ec ~ 0.2. For the state of two atoms in the first well, 
n = — 1, the single- atom (and two-atom) resonances occur 
at 5 = Ec/2 — 0.78E r (and S = 0) when the resonant 



p-4 



Many-body quantum coherence and interaction blockade 



conditions E*(n) = E*(n + 1) (and E*(n) = E*(n + 2)) 
are satisfied. In which, the recoil energy E r — h 2 / (2mX 2 ) 
with the mass m for a 87 Rb atom, the Planck's constant 
h, and the short-lattice wave-length A = 765 nm. 
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Fig. 4: Peaks of single-atom resonant tunneling and plateaus 
of interaction blockade, (a) The fluctuations of relative par- 
ticle numbers and (b) the expectation values of half relative 
particle numbers versus S/Ec- Solid, dash and dash-dot lines 
correspond to Ec/Ej = fOO, fO and f, respectively. 

The analysis presented above describes the system prop- 
erties at zero temperature. For finite-temperature sys- 
tems, assuming a canonical ensemble statistics, the equi- 
librium states are mixed states described by the den- 
sity matrix p(T) = exp(-/3H)/Tr[exp(-(3H)]. Here, 
(3 = l/(k B T), Tr[p(T)} = 1, k B is the Boltzmann con- 
stant, and H is the Hamiltonian matrix. The thermody- 
namic average and variance of an arbitrary operator Q 
are (Q) T = Tr[p{T)Q] and A(Q) T = Tr[p(T)(Q - (Q)) 2 ], 
respectively. 

For the systems of fixed parameters, we find that the 
coherence fluctuations grow with temperature. However, 
for the systems at fixed temperature, the parametric de- 
pendence of the coherence fluctuations relies on the val- 
ues of temperature. For low temperatures, the coherence 
fluctuations grow with Eq/Ej. This finite-temperature 
behavior is qualitatively consistent with the experimen- 
tal data on noise thermometry with two Josephson-linked 
condensates of Bose atoms [8]. For higher temperatures, 
the coherence fluctuations grow with Eq/Ej except for 
the anomalous behaviors near the value Ec/Ej ~ 4/iV 2 . 
For every eigenstate, we numerically find that the quan- 
tum fluctuations of coherence increase with Ec/Ej. This 
indicates that the monotonous growth of the coherence 
fluctuations vs Eq/Ej (except for the anomalous region) 
in finite-temperature systems is a signature of the quan- 
tum fluctuations. 

Using the mean-field Hamiltonian (2) with the same pa- 
rameters, we estimate the coherence fluctuations with the 
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Fig. 5: Quantum and thermal fluctuations of the phase co- 
herence vs. Ec/Ej for symmetric systems of TV = 100 and 
different values of temperature. Circles and dots correspond to 
the cases Ej = 1.0 and Ej — 0.5, respectively. Dashed lines 
are estimated with the classical Boltzmann theory. 

classical Boltzmann theory. For systems of low temper- 
atures and very small values of Eq/Ej, the classical es- 
timations show excellent consistence with the quantum 
ones. Otherwise, the classical results greatly differ from 
the quantum ones. This means that the mean-field Hamil- 
tonian works good for systems of low temperatures and 
very small values of Eq /Ej. Because the mean- field treat- 
ment uses the coherent state approximation and only the 
low-energy states for systems of small Ec /Ej are close to 
coherent states, the excellent consistence for systems of 
low temperature and small values of Ec / Ej is a natural 
result of the mean-field approximation. 

In Fig. 5, we show the quantum and thermal fluctu- 
ations of the phase coherence vs. Ec/Ej for symmet- 
ric systems of N — 100 and different values of T. In 
which, the dashed lines are obtained from the mean-field 
model (2) and the other data are calculated from the full 
quantum model (1). The region near Ec/Ej > 4/iV 2 
is the region of degenerated highest excited states which 
corresponds to the macroscopic quantum self-trapping in 
mean- field models [20,21,29]. The anomalous behavior 
of the coherence fluctuations in this region is caused by 
the occupation of the high excited states of two-fold de- 
generacy when the temperature becomes sufficiently high. 
We also calculate the averaged energy spacing for differ- 
ent parameters, and the results show a minimum in the 
crossover region between normal and self-trapping regions. 
For fixed and sufficiently high temperatures, the popula- 
tions of high-excited states in this crossover region will 
larger than those in other regions, and then the coherence 
fluctuations show an anomalous peak. 

In conclusion, we have studied the mechanisms of 
coherence fluctuations and interaction blockade in two 
Josephson-linked Bose-Einstein condensates with repul- 
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sive interaction. Our results for the coherence fluctuations 
and their resonant suppression are qualitatively consis- 
tent with several recent experimental studies of coherence 
fluctuations [6,8], conditional tunneling of single atoms 
and atom pairs [9], number squeezing [10], and squeezing- 
assisted long-time coherence [11]. We have revealed how 
the fluctuations depend on the system parameters and 
temperature. The different parametric dependencies of 
the coherence fluctuations and the number squeezing pro- 
vide a route for controlling the coherence time exper- 
imentally by adjusting the barrier height (junction en- 
ergy) [11]. In particular, we have shown that a compe- 
tition between the potential asymmetry and inter-particle 
interaction leads to the resonant tunneling and interaction 
blockade. The resonant effects explain consistently the 
recent observations of the conditional tunneling of single 
atoms and atom pairs [9]. More promisingly, the mecha- 
nisms of resonant single-atom tunneling and interaction 
blockade provide an alternative approach for designing 
single-atom devices for applications in sensitive metrol- 
ogy and information technology. At higher temperatures, 
an anomalous behavior of the coherence fluctuations is an 
indicative of the energy level degeneracy corresponding to 
the macroscopic quantum self-trapping effect [20,21]. 
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